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(Dated: February 7, 2008) 

Using the exclusive decay B 9 — * J/ip(fi + n~)<j>(K + K~), we report the most precise single 
measurement of the B 9 lifetime. The data sample corresponds to an integrated luminosity 
of approximately 220 pb _1 collected with the D0 detector at the Fermilab Tevatron Collider 
in 2002-2004. We reconstruct 337 signal candidates, from which we extract the B 9 lifetime, 
t(B 9 ) = 1-444+q; ,?, 8 (stat) ± 0.020 (sys) ps. We also report a measurement for the lifetime of the 
B° meson using the exclusive decay B° -» JW(/i + fi~)K* 9 (892)(K + ir~). We reconstruct 1370 
signal candidates, obtaining r(B 9 ) = 1.473l a050 (stat) ± 0.023 (sys) ps, and the ratio of lifetimes, 
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t(B°)/t(B°) = Q.980lg;gn ( stat ) ± °- 003 ( s y s )- 
PACS numbers: 14.40.Nd, 13.25.Hw 

Lifetime differences among hadrons containing b 
quarks can be used to probe decay mechanisms that go 
beyond the quark-spectator model Q. In the charm 
sector, lifetime differences are quite large 0; however, 
in the bottom sector, due to the larger 6-quark mass, 
these differences arc expected to be smaller. Phenomeno- 
logical models predict differences of about 5% between 
the lifetimes of B + and B°, but no more than 1% be- 
tween B° and B® lifetimes yj. These predictions are 
consistent with previous measurements of _B-meson life- 
times 0. It has also been postulated that the life- 
times of the two CP eigenstates (of the B®-B® system) 
differ. This could be observed as a difference in lifetime 
between B® semileptonic decays, which should have an 
equal mixture of the two CP eigenstates, and the lifetime 
for B® — > J/tpcp, which is expected to be dominated by 
the CP-even eigenstate 3] . 

In this Letter, we report a measurement of the life- 
time of the Bg meson using the exclusive decay channel 
B a s -> J/tp<p, followed by J/tp -> fj+ pT and (p -> K+K~. 
The lifetime is extracted using a simultaneous unbinned 
maximum likelihood fit to masses and proper decay 
lengths. We also measure the lifetime of the B° meson 
in the exclusive decay 1 B° -> J/tp if* (892), followed 
by J/tp -> fx + n~ and K*°(892) -> K + tt~ , and extract 
the ratio of the lifetimes of the B° a and B° mesons. The 
analysis is based on data collected with the D0 detec- 
tor in Run II of the Fermilab Tevatron Collider during 
the period September 2002-February 2004, which cor- 
responds to approximately 220 pb _1 of pp collisions at 
y/s = 1.96 TeV. 

The D0 detector is described in detail elsewhere 
We describe here only the detector components most rel- 
evant to this analysis. The central-tracking system con- 
sists of a silicon microstrip tracker (SMT) and a central 
fiber tracker (CFT), both located inside a 2 T supercon- 
ducting solenoidal magnet The tracking system and 
solenoid is surrounded by a liquid argon calorimeter. The 
SMT has « 800, 000 individual strips, with typical pitch 
of 50 — 80 /jm, and a design optimized for tracking and 
vertexing capability for |^| < 3, where 77 = — ln[tan(6>/2)] 
is the pseudorapidity and 9 is the polar angle measured 
relative to the proton beam direction. The system has a 
six-barrel longitudinal structure, each with a set of four 
layers arranged axially around the beam pipe, and inter- 
spersed with sixteen radial disks. The CFT has eight thin 
coaxial barrels, each supporting two doublets of overlap- 



Unless explicitly stated, the appearance of a specific charge state 
will also imply its charge conjugate throughout this Letter. 



ping scintillating fibers of 0.835 mm diameter, one dou- 
blet parallel to the beam axis, and the other alternating 
by ±3° relative to this axis. Light signals are transferred 
via clear light fibers to solid-state photon counters that 
have a quantum efficiency of approximately 80%. The 
muon system resides beyond the calorimeter, and consists 
of a layer of tracking detectors and scintillation trigger 
counters before 1.8 T toroidal magnets, followed by two 
similar layers after the toroids. Muon identification for 
1 77 1 < 1 relies on 10 cm wide drift tubes, while 1 cm wide 
mini-drift tubes are used for 1 < \rj\ < 2. Coverage for 
muons is partially compromised at the bottom of the de- 
tector where the calorimeter is supported mechanically 
from the ground. Luminosity is measured using plastic 
scintillator arrays located in front of the end calorimeter 
cryostat, covering 2.7 < \rj\ < 4.4. 

The data collection consists of a three-level trigger sys- 
tem, designed to accommodate the high luminosity of 
Run II. The first level uses information from the track- 
ing, calorimetry, and muon systems to reduce the rate 
for accepted events to « 1.5 kHz. At the next trigger 
level, with more refined information, the rate is reduced 
further to « 800 Hz. The third and final level of the 
trigger, with access to all of the event information, uses 
software algorithms and a computing farm and reduces 
the output rate to w 50 Hz, which is recorded for further 
analysis. We did not require the presence of any specific 
trigger in the event selection. 

Reconstruction of S° — > J/tp<p candidates requires 
a pair of oppositely charged muons that are identified 
by extrapolating charged tracks into the muon system 
and matching them with hits in the muon system. All 
charged tracks used in this analysis are required to have 
at least one hit in the SMT. We require that muon candi- 
dates each have a minimum transverse momentum px > 
1.5 GeV/c and that they form a common vertex, accord- 
ing to the algorithm described in Ref . |]| , which is based 
on a fit requiring a % 2 probability greater than 1%. The 
dimuon system was required to have an invariant mass 
between 2.90 and 3.15 GeV/c 2 and transverse momen- 
tum above 4.5 GeV/c. The dimuons are then combined 
with another pair of oppositely charged tracks, each with 
Pt > 0.8 GeV/c, consistent with the decay <j> — > K + K~ . 
The <p candidate was required to have an invariant mass 
between 1.008 and 1.032 GeV/c 2 and transverse momen- 
tum greater than 2 GeV/c. A four-track secondary ver- 
tex is fitted to the products of the J/tp and <p decays, 
and required to have a x 2 probability of at least 1%. 
The mass of the J/tp candidate is constrained in the fit 
to the world average J/tp mass of 3.097 GeV/c 2 @, the 
constraint does not take into account the uncertainty in 
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the J/tp mass. The resulting B® candidate is required to 
have pt > 6.5 GeV/c. We allow only one B° s candidate 
per event, and when multiple candidates exist, we choose 
the one with the best vertex probability. The resulting 
invariant mass distribution of the J/ijj-(j) system is shown 
inFig.OIa). 

Each primary vertex is reconstructed using tracks and 
the mean beam-spot position. The latter is determined 
for every data run, where a typical run lasts several hours. 
The initial primary vertex seed is constructed using all 
available tracks; a track is removed when it causes a 
change of more than 9 units in the \ 2 f° r a fit to a com- 
mon vertex. Thejjrocess is repeated until no more tracks 
can be removed 
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FIG. 1: (a) Mass distribution for B® candidate events. Points 
with error bars show the data, and the solid curve represents 
the result of the fit. The mass distribution for the signal 
is shown in gray; (b) same distribution after requiring the 
significance of the lifetime measurement to be ct/cj{ct) > 5. 

We take the four-track vertex as the position of the 
secondary vertex. To determine the distance traveled 
by each candidate, we calculate the signed transverse 
decay length (in a plane transverse to the direction of the 



beam), L x 



(Pt/pt), where x is the length vector 



pointing from the primary to the secondary vertex and pr 



is the reconstructed transverse momentum vector of the 
B®. The proper decay length of the B® candidate is then 
defined as cr = L xv {Mbo /Pt)> where Mgo is taken as 
the world average mass of the B® meson 5.3696 GeV/c 2 

a 

Figure ^b) shows the reconstructed invariant mass 
distribution of the B® candidates after a proper decay 
length significance requirement of ct/u{ct) > 5 is im- 
posed, where cr(cr) is the uncertainty on cr. The strong 
suppression of the background by this cut implies that 
the background is dominated by zero lifetime vertices, as 
expected. 

The proper decay length (without any restriction on 
significance) and the invariant mass distributions for can- 
didates passing the above criteria are fit simultaneously 
using an unbinncd maximum likelihood method. The 
likelihood function C is given by: 

N 

i 

where T s is the product of probability density functions 
for mass and proper decay length for Th is the equiv- 
alent for background, f s is the fraction of signal, and N 
is the total number of candidate events in the sample. 

The proper decay length for signal events is modeled 
by a normalized exponential-decay function convoluted 
with a Gaussian function of width equal to the uncer- 
tainty on the proper decay length, which is typically ~ 
25 jUm. This uncertainty is obtained from the full co- 
variance (error) matrix of tracks at the secondary vertex 
and the uncertainty in the position of the primary vertex. 
The uncertainty is multiplied by a scale factor that is a 
parameter in the fit to allow for a possible misestimate 
of the decay length uncertainty. The mass distribution 
of signal events is modeled by a Gaussian function. 

The proper decay length for the background is 
parametrized as a sum of a Gaussian function centered 
at zero and exponential decay functions, with two short- 
lived components and a long-lived term. The long-lived 
component accounts for heavy- flavor backgrounds, while 
the other terms account for resolution and prompt con- 
tributions to background. The mass distribution for the 
background is modeled by a first-order polynomial. 

To determine the background we use a wide mass range 
of 5.078-5.636 GeV/c 2 in the fit, corresponding to 4236 
B° candidates. The number of background candidates 
in this range is sufficiently large to measure the parame- 
ters of the background with high accuracy and therefore 
extract a good measurement of the signal fraction and 
ct(B®). The fit provides the cr and mass of the B®, the 
shapes of the proper decay length and mass distributions 
for the background, and the signal fraction. Table [i] lists 
the fit values of the parameters and their uncertainties. 
The distribution of proper decay length and fits to the 
B® candidates are shown in Fig. HJa) . 
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TABLE I: Values of the extracted mass Mb, resolution on the 
reconstructed mass gm, the measured ct, the signal fractions 
f s , and the scale factor s. 



Parameter 


B° a -» J/U 
fit values 


B° -> J/V'A'* (892) 
fit values 


Mb 
cm 

CT 

fs 
S 


5357.0 ± 2.5 MeV/c 2 
32.9^2.3 MeV/c 2 

433+2? V m 
0.0796 ± 0.0058 
1.142 ± 0.028 


5271.2 ± 1.5 MeV/c 2 
37.9±{\t MeV/c 2 

442±^ )im 
0.0446 ± 0.0018 
1.128 ± 0.009 




Proper Decay Length (cm) 




10 



-0.1 -0.05 



0.05 0.1 0.15 0.2 

Proper Decay Length (cm) 



FIG. 2: Proper decay length distributions for (a) B® and 
(b) B° candidates. The points with error bars show the data. 
The solid curve shows the total fit, the dashed curve the back- 
ground component, and the shaded region the signal. 



With a very similar four-track topology in the final 
state, the exclusive decay B° J /ipK*° (892) followed 
by J/ip -> fi+n- and if* (892) -> K+n- is recon- 
structed using the same selection criteria and algorithms 
as for the B® channel described above. The only dif- 
ferences are the requirement that the pr of the pion be 
greater than 0.5 GeV/c, and the selection of the if* (892) 
candidates. The combination of two oppositely charged 
tracks, assuming the pion mass for one and the kaon mass 



for the other, that gives an invariant mass closest to the 
mass of the iv*°(892) @ is selected for further study. 
The invariant mass of these combinations is required to 
be between 0.850 and 0.930 GeV/c 2 . Using the sample 
of B° candidates in the mass range 4.935-5.610 GeV/c 2 , 
corresponding to 30692 candidates, we determine the ct 
and mass of the B° using exactly the same procedure as 
used for B° s mesons. Results are also given in Table HJ 
and the distribution of proper decay length is shown in 
Fig. lib). 

Detailed Monte Carlo studies were performed on en- 
sembles of events comparable to data samples, with sim- 
ilar resolutions, pulls, fitting and selection criteria. No 
significant biases resulting from our analysis procedures 
were observed. To test the stability of the fit results for 
B® and B° mesons, we split each data sample into two 
roughly equal parts in order to study different kinematic 
and geometric parameters, compared the fit results, and 
found consistency within their uncertainties. We varied 
the selection criteria and mass ranges, and did not ob- 
serve any significant shifts. Using Monte Carlo samples 
with different input proper decay lengths in the range 
340 to 560 /mi, we checked the response of our fits to 
this variation, and found it to be linear in this range. We 
studied the contamination of our sample from cross-feed 
between B° and B° using Monte Carlo events. The es- 
timated contamination is 4.4% for B° s and 1.1% for B°, 
with invariant mass spread almost uniformly across the 
entire mass range. Therefore, their contributions are in- 
cluded in the long-lived heavy-flavor component of the 
background. To study possible biases from our fitting 
procedure, we used toy Monte Carlo ensembles with the 
same statistics as our data and with distributions match- 
ing those in data. These samples were fit, and the result- 
ing means and widths of the distributions of extracted 
parameters are consistent with the fits to data. 

Other sources of systematic uncertainty have been con- 
sidered, and the contributions are listed in Table ITT1 For 
the B° lifetime, there arc major contributions from de- 
termination of the background, the model for resolution, 
and the reconstruction of the secondary vertex. To deter- 
mine the systematic error due to the uncertainty in the 
background, we considered different models for the mass 
and decay-length distributions. In particular, to account 
for any model dependence on the invariant mass of misre- 
constructed heavy-flavor hadrons, we fit the probability 
distributions separately in the lower-mass and higher- 
mass side-band regions, and found the long-lived com- 
ponent to have different exponents. Combining the two 
lifetime values for the long-lived components, we modi- 
fied the functional form of the long-lived component for 
the global background in our fit. The two long-lived 
components were combined using a weighting parame- 
ter w = 0.98^Q'3g. This weighting parameter was varied 
by its uncertainty. The largest difference in the ct(B®) 
observed in these variations of background modeling was 
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found to be 4 /im, and is taken as the systematic uncer- 
tainty due to this source. The effect of uncertainty in the 
proper decay length resolution was studied by using an 
alternative resolution function consisting of two Gaussian 
functions (with the same mean but different width), re- 
sulting in a difference in the fitted ct(B®) of 3 /im. Uncer- 
tainty or biases in the determination of the secondary ver- 
tex were estimated using secondary vertices constructed 
with the J/tp tracks only, resulting in a ct(B®) shift of 
3 /im. The contribution from the uncertainty on the de- 
tector alignment is estimated by reconstructing the B® 
candidate events with the position of the SMT sensors 
shifted radially outwards by the alignment error in the 
radial position of the sensors. The resulting difference in 
fitted proper decay length of 2 /im is taken as the sys- 
tematic uncertainty due to possible misalignment. The 
total systematic uncertainty from all these sources added 
in quadrature is 6 /im. The systematic uncertainties in 
the measurement of the ct(B°) are determined in the 
same way as for the B®, and each contribution is listed 
in Table HD 



TABLE II: Summary of systematic uncertainties. 





ct(B°) 

H 


ct{B°) 
(/im) 


r(B° s )/r(B°) 


Alignment 


2 


2 


0.000 


J/ip vertex 


3 


4 


0.002 


Model for resolution 


3 


3 


0.000 


Background 


4 


5 


0.002 


Total 


6 


7 


0.003 



We determine the lifetimes of the B® and B° mesons, 
t(B q s ) = 1.444tS ( stat ) ± °- 020 (sys) PS, 



uncertainty of that source, as shown in Table HU The 
sum in quadrature of all contributions is reported as the 
overall systematic uncertainty on the ratio of lifetimes 
including correlations between the two lifetime measure- 
ments. 

In conclusion, we have measured the B° s and B° life- 
times in exclusive decay modes in pp collisions. The mea- 
surements are consistent with previous results The 
value of the B® lifetime obtained in this analysis is the 
most precise measurement from any single experiment. 
The ratio of the lifetimes is also in good agreement with 
QCD models based on a heavy quark expansion, which 
predict a difference between B® and B° lifetimes of the 
order of 1% Q. 
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t(B°) = 1.4731^50 (stat) ±0.023 (sys) ps. 

Both results are consistent with the current world 
averages of t(B2) = 1.461 ± 0.057 ps and t(B°) = 
1.536±0.014 ps |2|. We note that measurements using B® 
semileptonic events, where there is an equal mixture of 
CP-even and CP-odd states, dominate the current world 
average, while B® — > J/ip<fi has a different composition of 
CP-even and CP-odd states as discussed earlier 

Using our results we determine the ratio of B® / B° life- 
times to be 

^0- = 0.980 ±°;°™ (stat) ± 0.003 (sys), 

where statistical uncertainties were propagated in 
quadrature, and the systematic uncertainty was evalu- 
ated by adding each contribution to the corresponding 
central value, and evaluating a new ratio, with the dif- 
ference from the nominal value taken as the systematic 



[*] Visitor from University of Zurich, Zurich, Switzerland, 
[t] Visitor from Institute of Nuclear Physics, Krakow, Poland. 

[1] M. B. Voloshin and M. A. Shifman, Sov. Phys. JETP 
64, 698 (1986); I. Bigi and N. G. Uraltsev, Phys. Lett. 
B 280, 271 (1992); I. Bigi, Nuovo Cimento A 109, 713 
(1996); F. Gabbiani, A. I. Onishchenko and A. A. Petrov, 
arXiv:hep-ph/0407004 

[2] S. Eidelman et al. (Particle Data Group), Phys. Lett. B 
592, 1 (2004). 

[3] I. Bigi et al. in "B Decays" , 2nd edition, edited by S. Stone, 
(World Scientific, Singapore 1994) p. 13 ; I. Dunietz, Phys. 
Rev. D 52, 3048 (1995); M. Beneke, G. Buchalla, I. Duni- 
etz, Phys. Rev. D 54, 4419 (1996). 

[4] T. LeCompte and H. T. Diehl, "The CDF and D0 Up- 
grades for Run II", Ann. Rev. Nucl. Part. Sci. 50, 71 
(2000); V. Abazov, et al, in preparation for submission 
to Nucl. Instrum. Methods Phys. Res. A. 

[5] J. Abdallah et al. (DELPHI Collaboration), Eur. Phys. J. 
C 32, 185 (2004). 

[6] The exact fraction of CP-even and CP-odd decays is un- 



8 



known. In this analysis, from Monte Carlo, the relative 
efficiency for CP-even and CP-odd states is 0.99 ± 0.01. 



